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In critical phenomena, singular behaviors arise not only for thermodynamic quantities but also 
for transport coefficients. We study this dynamic critical phenomenon in the AdS/CFT duality. 
We consider black holes with a single R-charge in various dimensions and compute the R-charge 
diffusion in the linear perturbations. In this case, the black holes belong to model B according to 
the classification of Hohenberg and Halperin. 
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I. INTRODUCTION 

The AdS/CFT duality is a powerful tool to study gauge theories at strong coupling. At zero 
temperature, there are many circumstantial evidences for the duality, but the finite temperature, 
dynamic cases are less understood. In fact, the results for these cases such as the universality of fj/s 
(fj: shear viscosity, 1 s: entropy density) are mostly regarded as predictions for dual gauge theories. 
(See Ref. [l| for a review.) Obviously, the problem comes from the fact that it is extremely hard to 
compute gauge theories at strong coupling. 

The critical phenomena may be useful to check some aspects of the duality. At the critical point, the 
correlation length £ diverges, and only a few low-energy, long- wavelength modes become dominant. 
As the consequence, it is possible to capture the behavior of physics which does not depend on the 
details of microscopic physics. This property enables one to check the AdS/CFT duality even if one 
cannot compute gauge theories at strong coupling. 

At the equilibrium, the effect of diverging correlation length appears as the divergence of thermody- 
namic quantities such as the specific heat (static critical phenomena). The divergence is parametrized 
by the static critical exponents. For normal statistical systems, the effect of diverging correlation 
length also appears in nonequilibrium cases (dynamic critical phenomena). For example, the relax- 
ation time 2 and the other transport coefficients diverge, which is known as the critical slowing down. 
The divergence is parametrized by the dynamic critical exponents, and these dynamic critical expo- 
nents are related to the static critical exponents. Our aim is to see these singular behaviors from 
black hole physics and to see if they agree with the predictions from the theory of dynamic critical 
phenomena. 

The simplest AdS/CFT duality involves the Schwarzschild-AdSs (SAdSs) black hole (with planar 
horizon), which has no phase transition at finite temperature. The SAdSs black hole with compact 
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1 In this paper, we use fj for the shear viscosity to avoid confusion with a static critical exponent n (See Sec. Ill Al . 

2 One should not confuse this relaxation time with the one appeared in the second-order hydrodynamics 
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horizon such as S 3 can have a phase transition, but it is not suitable for our purpose: First, it is a 
first-order phase transition; Second, the theory has no hydrodynamic limit [ui(q) —> as q — > 0] since 
q takes only the discrete values. Similarly, the Reissner-Nordstrom-AdSs black hole with compact 
horizon has a critical point similar to QCD, but it is not suitable since it has no hydrodynamic limit. 
A simple system which has both of the second-order phase transition and the hydrodynamic limit is 
the D = 5 single R-charge black hole with planar horizon 0,0)13]. 3 

In this paper, we consider single R-charge black holes in various dimensions. We consider the linear 
perturbations in the backgrounds, compute the relaxation time for the R-charge diffusion, and show the 
critical slowing down. We show that the black holes belong to model B according to the classification 
of Hohenberg and Halperin Q, i.e., the same universality class as the uniaxial ferromagnet. (One can 
regard our systems as model H, the same universality class as the liquid-gas phase transition. Sec 
Sec. EH) 

In the next section, we briefly review both static and dynamic critical phenomena [7[ . We analyze 
the dynamic critical phenomena for the D = 5 single R-charge black hole in Sec. IIIII In Sec. IIVI we 
extend the analysis to the D = 4, 7 single R-charge black holes. We discuss some related issues in 
Sec. El 4 

II. CRITICAL PHENOMENA 
A. Static critical phenomena 

For the second-order phase transition, thermodynamic quantities diverge as powers of the correlation 
length £, and these powers, static critical exponents, are universal. Namely, they depend on the 
symmetry, spatial dimensionality, and so on, but not the other details of the interactions. Different 
physical systems may belong to the same universality class; e.g., the 3d Ising model and the liquid-gas 
phase transition. These critical exponents also satisfy static scaling relations, which implies that not 
all critical exponents are independent. In fact, these scaling relations are derived from a scaling law 
for the thermodynamic potential. 

As an example, consider the ferromagnetic phase transition. In this case, the magnetization m and 
the external magnetic field h are the order parameter and a control parameter, respectively, which are 
related by 




(11.1) 



where f2 = £l(T, h) is the Gibbs free energy whose variation is given by 

fl = -sdT-mdh . (11.2) 



3 In this paper, we denote D as the number of bulk (noncompact) spacetime dimensions and d as the number of 
boundary spatial dimensions. 

4 Recently, the dynamic critical exp onent has been discussed in the context of the Schrodinger group in Refs. |25l l26t 
I27ll28t l29l. l3(](l . See also Refs. [3l[. 1 32il which overlap with these works. 
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The static critical exponents (a, /?, 7, S, v, rj) are denned by 

fd 2 fl\ 

specific heat: C h := -T — — oc \e T \ 01 , (II.3a 



\dT 2 J 



h 



spontaneous magnetization: to oc \ct\ {T < ^c), (II. 3b) 

/ ^777, \ 

magnetic susceptibility: %t := ygfo ) K l £T l~ 7 > (II. 3c) 

critical isotherm: to oc |/i| 1/l5 (T = T C ), (II.3d) 

correlation function: G(r ) oc e~ r/ * (T^T C ), (II.3e) 

ocr- d+2 -" (T = T C ), (H.3f) 

correlation length: £ oc |£t| _1/ , (H-3g) 

where £t := (7 1 — T c )/T c , and 0? denotes the number of spatial dimensions. The correlation function 
is dchned by 

G(r) := (to(F)to(O)) oc . (IL4) 

By dehnition, \t oc G(q = 0), where G(q) is the Fourier component of G(r). Since G(g) oc q~ 2+v oc 

£ 2_r? (?<z) _2+,7 j the susceptibility diverges at the critical point since the correlation length diverges. In 
fact, the scaling law below claims that the divergence of all thermodynamic quantities is due to the 
diverging £. 

There are 6 static critical exponents in Eq. (|II.3|) . but not all are independent, and they satisfy 
static scaling relations: 

a + 2/3 + 7 = 2, (II.5a) 

7 = (3(6 - 1) , (H.5b) 

7 = 1/(2 - 77) , (II.5c) 

2 - a = . (II.5d) 

A relation involving d such as Eq. (|II.5dp is known as a hyperscaling relation. Because of these 
relations, only 2 are independent among 6 exponents, which suggests that there is some structure 
behind these relations. 

In fact, these relations can be derived from the static scaling law for the free energy: 

fi(e T) h) = b- d n(b yT e T , b Vh h) . (II.6) 

Equation (|II.6[) implies 

n(e T , h) = e d T /vT n(l,e T Vh/yT h) =: 4 /yT V{e T Vh/vT h) . {11.7) 
The scaling law gives a scaling law for the correlation function: 

G{q, e T , h) = b 2yh - d G(bq, b yT e T , b Vh h) . (II.8) 
The static scaling law determines 6 critical exponents in terms of 2 parameters (yr,y?i): 

a = 2 - d/y T , (II.9a) 

(3 = , (II.9b) 

VT 

7 = , (IL9C) 
Vt 

5 = , (II.9d) 
a- Vh 

v = l/VT , (H.9e) 

Tj = d-2y h + 2 . (II. 9f) 
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Equations ()II.9a|) - (|II.9fj) satisfy the scaling relations (|II.5a[) - (III.5dp . It is useful to rewrite the scaling 
law (|IL8|) in a form similar to Eq. (|II.7[) : 

G(q,z,h)=e- n $m,e h h) , (mo) 

where Eqs. (|II.9el) and (|II.9fj) arc used. Note q appears in the form (£<t): this is a typical form in a 
scaling law. 

As an example, consider the (uniaxial) ferromagnets such as the Ising model. In the mean-field 
approximation, the Ising model is described by the Ginzburg-Landau theory: 



/[to; T,h] = J dx 



c . , n Oi 9 b a 

— ( vm) + — to + — to H — mn 

2 2 4 



(11.11) 



where a = a (T — T c ) + • • • , b = b + • • • , and c = c + • • • (The dots represent the terms which 
do not contribute near the critical point). As usual for the Legendre transformation, I[m;T, h] is a 
pseudofree energy, and the free energy is given by eliminating to: Q(T,h) — min m I[m; T, h]. The 
static critical exponents for the Ginzburg-Landau theory are 

{a, 0,"f, 5, is, rj) = k^l,3,i 0J (11.12) 

for d = 3. These exponents satisfy the scaling relations except the hyperscaling relation (|II.5d[) . In 
fact, it is easy to show that the pseudofree energy pi. lip satisfies the scaling law ()II.6[) . 
There are some caveats which are relevant to our later discussion: 

• The mean- field theory ignores the effects of statistical fluctuations. The Ginzburg criterion tells 
that the mean- field theory is trustable for d > d c , where d c is the upper critical dimension given 

by 

dc = 2 i±^ = az£ . (II .i3) 

V V 

Namely, the static critical exponents obtained from a mean-field theory is unreliable below d c . 
For the Ginzburg-Landau theory, d c = 4. On the other hand, we consider the large-iV c limit, 
where the effects of fluctuations are suppressed. 

• The hyperscaling relation (|II.5d[) may fail above d c (or even below it in some problems). This is 
because the scaling law for the free energy pi.6[) is modified by the so-called dangerous irrelevant 
operators. 



B. Dynamic critical phenomena 



In the dynamic case, the order parameter varies slowly in time. Near the critical point, a large 
domain is favorable due to the diverging correlation length. Because of the large size, it takes a long 
time for thermal fluctuations to flip spins for a whole domain. This is known as the critical slowing 
down, and the divergence of the relaxation time is parametrized by the dynamic critical exponent z. 

The idea of scaling law can be naturally extended into the dynamic case. The static scaling law for 
the correlation function (|II.10[) is extended into the dynamic scaling law: 

G(u,q,e T ,h) = b 2 - r ' +z G(b z uj,bq,b yT e T ,b Vh h) (11.14) 
= ?-v+>S>(Z' U) tf,?»h) , (11.15) 

where z is a dynamic critical exponent. The scaling law suggests that the relaxation time diverges 
near the critical point: 

r q ^ Epical ~ ^ ^ (T^T C ), (II. 16a) 

~q- z {iq) z ^q) <xq- z (T = T C ). (II.16b) 
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[In reality, the fluctuation-dissipation theorem is necessary to show Eqs. (|H.16[) . and we skip some 
steps here.] 

The dynamic universality classes were classified by Hohenberg and Halperin [|| : they are known as 
model A, B, C, H, F, G, and J. These models are further classified by the values of dynamic critical 
exponents. The dynamic universality class depends on additional properties of the system which 
do not affect the static universality class. In particular, conservation laws play an important role 
to determine dynamic critical exponents. A conservation law forces the relaxation to proceed more 
slowly. As a consequence, even if two systems belong to the same static universality class, they may 
not belong to the same dynamic universality class. 

We consider model B example below. (See Sec. IV Al for model H.) Model B corresponds to a system 
with a conserved charge. This is the case for uniaxial ferromagnets, but we use the conventions for the 
charge diffusion below for a later purpose. In the hydrodynamic limit, it is governed by the diffusion 
equation: 

= d tP (t,x)-D R Ap(t,x) ^ 0=(u + iD R q 2 ) p(u,q) , (11.17) 

where D R is the diffusion constant. Note that the diffusion constant D R is related to the susceptibility 
Xt and the conductivity A as 



X = D RXT , (11.18) 



from the definitions of A and D R : 



J, = — AV^/j = -A I ^ ) Vip , (11.19) 



dp , 

J, = -D R W iP . (11.20) 

It is important to distinguish three different regimes: the hydrodynamic regime (£g <C 1), the 
critical regime (£q 3> 1), and the overlapped regime (£g ~ 1). 

1. The hydrodynamic regime: In this regime, the charge diffusion is governed by the diffusion 
equation, and the relaxation time T q for wave number q is given by 

l/r q ~ -9f(w) = D R q 2 . (11.21) 

For model B, A is nonsingular on the critical point. Because xt diverges on the critical point, 
this means that D R vanishes on the critical point, namely the critical slowing down. Then, 
Eqs. (|HTT8|) and (|Tl3T|) give 

1/ Tq oc xt 1 I 2 • (H.22) 

2. The critical regime: On the other hand, in the critical regime, higher powers of q in the dispersion 
relation are no longer negligible. In this regime, the dynamic scaling relation (|II.16|) dictates the 
dispersion relation l/r q oc q z . 

3. The overlapped regime: These two behaviors should match smoothly at ^gft! 1. Thus, Eq. pi.22[) 
should be rewritten as 

l/r q cx C (i - V) m 2 , (H.23) 

where a static scaling relation is used. 

Comparing Eq. pi.23|) with Eq. (|II.16a|) . one gets 

z = 4 - 77 . (11.24) 

Note that Eq. (|II.24[) relates the dynamic critical exponent z to the static critical exponent 77. 

The important assumption of model B is that A remains finite on the critical point. We will check 
that this is indeed the case for single R-charge black holes. 
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III. DYNAMIC CRITICAL PHENOMENA FOR THE D = 5 SINGLE R-CHARGE BLACK 

HOLE 

In this section, we show that the critical slowing down indeed occurs at the second-order phase 
transition for the D = 5 single R-charge black hole [3| . 



A. The D — 5 single R-charge black hole and its thermodynamics 

In the simplest AdS/CFT duality, the gravity dual is SAdSs x S 5 . One can add an angular momen- 
tum along S 5 , which is known as the "spinning" D3-brane solutions [J, [f| . The angular momentum 
becomes a charge after the S 5 reduction. The symmetry of S 5 corresponds to an internal symmetry of 
the Af — 4 super- Yang-Mills theory, R-symmetry SU(4:)r. The R-symmetry group SU(4)r is rank 3, 
so one can add at most three independent charges. The three-charge solution is known as the STU so- 
lution [H . When all charges are equal, the STU solution is nothing but the Reissner-Nordstrom- AdSs 
black hole. 

The S 5 reduction of D — 10 type IIB supergravity gives D = 5 Af = 2 gauged U(l) 3 supergravity 
0: 



r-\ Ga f; u F j V-\ Gij (v M x<) (v^) 



4 

1? 



3 1 
Y — 

1=1 



1 

24 



e - u F l F j A\ 



(III.1) 



Here, we use the unit 2 k| = l&itGs = 1, and are the field strength of the 3 U(l) gauge fields 
A^ (i = 1, 2, 3). The fields X 1 (i — 1, 2, 3) represent 3 real scalar fields which are not independent but 
are subject to the constraint X 1 X 2 X 3 = 1, and their moduli space metric Gij is given by 



G 



diag[ {X'Y\ {X*)-\ (X 3 ) 



^3\-2 



For the single R-charge black hole, use the ansatz 

F% = f% = , X 2 = X 3 = 1/Vxi =: H 1 ' 3 

which is consistent with the equations of motion. Then, the action becomes 

"2 1 i\irj\2 



C 5 



-9 



8 3 H 2 



2V 



where 



2V = ( H 2/3 + 2H~ 1/3 



A^ := (V2/i) A 1 ^, and F^ := 2d\ LlM A v y The equations of motion are given by 

iV = 3 (V ^^ ( 2 V - g) + ^ H-' 3 F^F* - 9 -f{j H^ 3 F 2 + 4 
V„ ( H A ' 3 F^ v \ = , 



UY*H = -ZH— + — H A ' 3 F 2 
dH 4 



(III.2) 
(III.3) 

(III.4) 

(in.5) 

(m.6a) 
(III.6b) 
(III.6c) 
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and the single R-charge black hole (with planar horizon) is given by 

dsj = {nT ° L)2 H 1 ' 3 (- L dt 2 + dx 2 + dy 2 + dz 2 )+ -4^ # 1/3 du 2 , (III.7a) 



A„ = 7T To V2k(1 + k) jj (dt)„ , (III.7b) 

H = l + «u, (III.7c) 

/ = H - (1 + k) u 2 = (1 - u) {1 + (1 + k) u} , (HI.7d) 

where u £ [0, 1], with u = 1 corresponding to the location of the horizon. 

The black hole is characterized by two parameters (T , k) : the former gives the Hawking temperature 
for the neutral black hole, and the latter gives the R-charge. Thermodynamic quantities for the black 
hole are given by 

o 2 N 2 T 4 2 N 2 T 4 

e= N 8 cT ° + P= L^ASL {1 + K)j (111.8a) 

s=^^-VT+H, T=-^±^T , (III.8b) 

I 2 v 1 + k 



7r AT 2 T 3 , 19k 

p = ^ jl a/2«(i + k) , /1 = 7rTo ViT^' (m ' 8c) 

which satisfy the thermodynamic relations: 

de = Tds + pdp , dP = sdT + pdp , e + P = sT + pp. (III. 9) 

Hereafter, we study the black hole as the grand canonical ensemble, i.e., the system is characterized 
by (T,/x) not by (T ,k). The appropriate thermodynamic potential is the Gibbs free energy il — 
—P = e — Ts — pp. The quantities T and p satisfy 

2 7T T I K I 2 

= x - + J->2 (The equality holds for k = 2) . (III. 10) 

p y 2 y k 

For a given pair of (T, //) , there are two black hole solutions. These solutions are related by the 
transformation 



— t °\/I(TT^- (m - n) 

(This transformation is valid for D = 4, 7 single R-charge black holes in Sec. lIVI as well.) The solution 
thermodynamically realized is the one with the lower Gibbs free energy, which is the branch of k < 2. 
Thus, we focus on k < 2. 

Thermodynamic stability is determined from the behavior of thermodynamic fluctuations, e.g., the 
R-charge susceptibility and the specific heat at constant p: 

„^(V\ = m.i±pz±, (m . 12) 



dp 



T 



r-r( ds \ K 2.y2 T 3 (3-«)(2 + «) 

^-r^-j^-. N C T VT _ {2 _-, (111.13) 

which diverge at K = 2. On the other hand, the first derivatives of the Gibbs free energy such as s and 
p are regular there, so this is a second-order phase transition. Note that k = 2 imposes one condition 
in the (T, /x)-phase diagram, so n — 2 defines a critical line instead of a critical point. 
The thermodynamic quantities ()III.8[) determine the static critical exponents [H, Q : 

(a,M,5,v,ri)= (1,1 1,2 1 l) . (111.14) 



.2' 2' 2' '2 

(See App.H) 
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B. Critical slowing down for the D — 5 single R-charge black hole 

The single R-charge black hole has a second-order phase transition at k — 2 as we saw in the last 
subsection, so one expects a critical slowing down, where a relaxation time diverges. According to 
the discussion in Sec. Ill B[ if the single R-charge black hole belongs to model B, A remains finite on 
the critical line. In this subsection, we consider the linear perturbations for the black hole, and we 
explicitly show that A remains finite on the critical line. 

The R-charge conductivity A can be evaluated from a Green-Kubo formula: 



G i J x ) (cj,q = 0) 



A = — lim 

= Urn / dte 1 ^ fdx([J x (t,x), .4(0,0)]) , (111.15) 
2uj J_ OQ J 

where G x ^(u),q) is the retarded correlator of the R-charge current: 

G x f(uj,q) := -i f dt e iwt 0(t) f dx e~^ s {[J x (t,x), 4(0,0)]) . (HI.16) 



We evaluate this correlator from the AdS/CFT duality. 5 Since Eq. (|III.16|) is a vector mode correlator, 
it is enough to consider the vector modes of the bulk fields. 6 The perturbations from the background 
fields g , Afj, are denoted as 

h^u ■■= g^u - g^ v , a M := Ap- . (111.17) 

Since we consider the vector modes, the scalar field H does not fluctuate, i.e., H = H . Thus, we 
consider the perturbations which take the form 

h tx =: g xx T{u) e -*"*-H** , h zx =: g xx Z(u) e -*"*+** , a x =: | A(u) e ~^ t+ ^ , (111.18) 
with the other h^u = 0. 

Moreover, we need only the q = correlator since we evaluate the Green-Kubo formula (|III.15|) . 
This gives a further simplification, namely it is enough to consider the q = limit of the equations of 
motion for (T, Z,A): 

Z" + ^-lz' + ^-Z = , (III.19a) 

uf uf 2 

T ' + S A = °' (IIL19b) 

A " + (tf + j) A ' + A + 2 jw r = ' (IIL19c) 



5 Our results in this section can be obtained from the results of Ref. [To| (See Sec. 4 of the paper. See also Refs. [HE1 
I l.'l for related issues), but the computation presented here is easier than their computation. Both our computation 
and their computation consider the vector modes. They have done the computation to obtain the shear viscosity fj 
whereas we are interested in Dr. In the vector modes, the quasinormal mode method is used to obtain fj whereas 
the Green-Kubo formula method is used to obtain Dr. As the consequence, the former needs the 5^0 correlator, 
and the latter needs only the q = correlator. Our result is a limit of their result, but the computation is far easier. 

6 The electromagnetic perturbation or the R-charge current is decomposed as the vector mode and the scalar mode. 
Similarly, the gravitational perturbation or the energy-momentum tensor is decomposed as the tensor mode, the 
vector mode, and the scalar mode. 
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where 

W ■■= 2k • ( IIL2 °) 

In the q = limit, the field Z decouples from A and T, and Eqs. (|III.19b[) and piI.19c|) give an 
equation for A only: 



A "+\-n- + i j) A'+[—-^—t 1 u)A = Q. (111.21) 



Solving this equation is essentially enough to get the R-charge current correlator. 

The correlator is obtained from the Lorentzian prescription of the AdS/CFT duality [l4T ]: 

• Step 1 : Solve the bulk equations of motion by imposing the "incoming wave" boundary condition 
at the horizon and by imposing the boundary values of the bulk fields. The incoming wave 
boundary condition corresponds to the retarded condition, and the boundary values play the 
role of source terms for the boundary fields in the dual gauge theory. 

• Step 2: The total action consists of the bulk action pil.4[) . the Gibbons-Hawking term, and 
conterterm action which cancels divergences. Substitute the solution obtained in Step 1 into the 
action, and evaluate the boundary action quadratic in the boundary values. 

• Step 3: The correlator is obtained from the second derivative of the boundary action with respect 
to the boundary values. 

The first is to solve the bulk equation of motion (|III.21|) . which can be done perturbatively in to. 
Incorporating the incoming wave boundary condition at the horizon, we get 

o.(«) = «£»> £p {i + ^ + x {^h + 2 ~? (2+Ku)Hu+H) ) } + 0( " 2) ' (IIL22) 

where a x °^ represents the boundary value of a x (u): 

lim a x {u) = 4 0) . (111.23) 

The second is to obtain the boundary action. The boundary action is given by 

N?T 2 



N 2 T 2 f 

^boundary = ' lim / dt d % (/ H d x a x H ) 

lb u— >oj 



where the dots represent the terms which are not quadratic in a§ ■ Finally, following Step 3, we get 

16 64 7TV1 + K 



A = — lim 



G i * ) (L;,q = 0) 



(111.26) 



(k + 2) 2 N*T 

~o u 64 7rVl + k ' 

so A is finite on the critical line n — 2 as promised. 

Using Eq. (|III. 12|) . the R-charge diffusion constant is given by 

1 (I + k/2) 3 2-k ,„ 
D R = 1 1 ' ^ , (111.27) 

so Dr = on the critical line. 

To summarize, the dynamic universality class for the single R-charge black hole is model B with 
dynamic critical exponent z = 4 — 77 as predicted from the theory of dynamic critical phenomena. 
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IV. SINGLE R-CHARGE BLACK HOLES IN OTHER DIMENSIONS 



The extension to other dimensions is straightforward and similar. We consider the D = 4, 7 cases, 
which are obtained from rotating M2 and M5-branes, respectively. For the actions and the solutions, 
see, e.g., Ref. [f|. 

A. D — 4 single R-charge black hole 



For a single R-charge black hole, the D = 4 action is given by 

'2 q (T7U\2 



-4 



R _V HV2F2 _HVHf 

8 8 H 2 



2V 



(IV.l) 



where 



The single R-charge solution is 



2V = ^{H 1 ' 2 + H- 1 ' 2 ) 



9tr \ H 



Aft = -7tT V2k(1 + k) — (dt)/. , 

o rz 

= 1 + K U , 

/ = (I - u) {l + (I + k)u + {I + k)u 2 } 



which satisfies the equations of motion from Eq. (jIV.ip . 
Thermodynamic quantities are given by 



V2n 2 



Nl' 2 T* (1 + k) 



s = V2ir 2 (jj Nl /2 T 2 VTTk , 
p = V2n Nl' 2 T 2 ^2k{1 + k) 



P = 
T = 
/' 



V2 7T 2 

3 

\/l + K 



fit 



■H 1 ' 2 du 2 



Nl' 2 T* + 



2k 
1 + T 



4nT 2 k 



3 



1 + K 



from which one gets 



(IV.2) 



(IV.3a) 

(IV.3b) 

(IV.3c) 
(IV.3d) 



(IV.4a) 
(IV.4b) 
(IV.4c) 



ttT _ V3 



2 \ > V3 
3k I ~ 2 



(The equality holds for k = 3/2) 



(IV.5) 



Again there are two solutions of k for a single pair of (T, //), and we focus on k < 3/2. The R-charge 
susceptibility and the specific heat at constant /i are given by 



XT 



Nc /2 T 1 + 2k 
6y/2 3 -2k ' 



C, = V2n 2 I) NV*T* 



(2-k)(3 + 2k) 
VTTk"(3-2k) 



(IV.6) 
(IV.7) 



Thus, k = 3/2 is the critical line. The static critical exponents (a,/3,j,6) are the same as the D = 5 
case, and (u, rj) = (3/4, 4/3) from Eq. (|XTT|) @, 0. 
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Similar to the D = 5 case, the equations of motion for (T, Z, A) in the q — limit are given by 

(IV.; 



„ uf'-2f . 9w 2 H 
Z" + — — ^- Z' + — — Z = , 



uf 

KU 2 . 

T+ 2H A =° 



4/2 



A" 



A , i ^ 2h a 1 2(i+.) 

H + f ) A + A p A + fH -°" 



Equations (|IV.8|) give an equation for A only: 



A" 



" / 



A' 



9»»H *(! + *) ^ 



4/2 



/if 



Then, the solution is given by 



f-imT /2T 
^ =:a (0)Z_^ A + O(co 2 ) 



H 



A = l 



2ku 



2k 3 u 



3 + 2k u 



3(3-K)yTT^ 4(3-k) 2 (3 + 2k) 



3V^TT(3 - k) 2 In 



1-lt 



2^3 - k{k(k + 12) + 9} arctan 



y/(3 -«)(! + «) 
2 + (1 + k)u 



The boundary action is given by 



'boundary 



7Vc /2 Tq 

36%/^ 



lim / dt dx if H a'a x 



(IV.8b) 
(IV.8c) 

(IV.9) 
(IV.10) 



(IV.ll) 



(IV.12) 



Substituting Eq. pV.lip into the boundary action, we get the R-charge current correlator G xx (k>, q 
0), the R-charge conductivity A, and the R-charge diffusion constant Dr: 



GW(w,q = 0) = 



k Nc /2 T (3 + 2k) 2 iV c 3/2 



54^/2 



6 3 7rv/2(l + K) 



A 



Dt 



(3 + 2k) 2 N, 



■ \2 AT 3 / 2 



6 3 tt\/2(1 + k) 

3 (1+2k/3) 3 (1-2k/3) 



4ttT (1 + k)(1 + 2k) 
Thus, A is finite and Dr = on the critical line K = 3/2. 



(IV. 13) 
(IV. 14) 
(IV.15) 



B. D = 7 single R-charge black hole 

For a single R-charge black hole, the D = 7 action is given by 

3 (VJT) 2 



£7 = V~9 



R -^H e / 5 F 2 - 

8 10 H 2 



2V 



(IV.16) 
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2V = ± #4/5 + 6ff -i/5 ) _ (IV17) 



where 

9V — , _ 

I? I 2 

The single R-charge solution is 

\2 / f \ jl 



dg2 = 4(7tT L) gl/8 /_ /_ ^ 2 ^ ^ rfz 2 \ + 1/5 d „2 (I y 18a) 

' 9u V H J 4/u 2 



2 



A„ = -7rT 0V /2 K (l + «) — (dt)„ , (IV.18b) 

o -H 

i/ = i + ku 2 , (rv.isc) 

/ = (1-u){1 + m+(1 + k)u 2 } , (IV.18d) 

which satisfies the equations of motion from Eq. (lIV.16j) . 
Thermodynamic quantities are given by 

e=^( 2 S N?T$(1 + K ), P=^ fcY N*1${1 + k), (IV.19a) 



3tt 3 (0 iV c 3 T 5 VTT^, 



T ^7fe( 1 + i)' (IVJ9b) 



p = tt 2 ( | ) iV c 3 T 5 V2-(l + -) , M = —7r- \h^~ , (TV.19C) 



6 

3 J ~ c -u v---v- • , r- 3 y 1 + K 

from which one gets 

if = ^ vf (Vf + \/f ) - vf (The equaiity hoids f ° r k = 3) • (iv ' 2o) 

Again there are two solutions of K for a single pair of (T,fi), and we focus on k < 3. The R-charge 
susceptibility and the specific heat at constant /i are given by 

c «'^(S) NiT °wm0k- ,,v ' 22) 

Thus, k = 3 is the critical line. The static critical exponents (a,/3, 7, S) are the same as the D = 5 
case, and (v,rj) = (3/10, 1/3) from Eq. (fATTTI) 

The equations of motion for (T, Z, A) in the q = limit are given by 

Z» + Zf-^* + ¥4z = 0, (IV.23a) 

2 

f + ^-A = 0, (IV.23b) 
H 



... 1 H' f 1\ . 9to 2 H . Auil + n)^. „ 

A + -EF + - - A + ~r^T A + T = • IV.23c 



Equations (|IV.23|) give an equation for A only: 



,„ , H' f 1\ „. /9to 2 H 4k(1 + k) , , 
1 ( H + 7 " «) ^ + l^F /I/ 2 ^ 1 1 : : H ' 
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Then, the solution is given by 

A=:a^ J —— A + 0(u; 2 ) , 



H 

KU 2 

A = 1 + — 



(IV.25) 



6(1 + re) (3 + ku 2 ){2(k - 3)re - 9) arctan- 



12VTT^(3 + k)(3 + 4k) 3 / 2 y v 2 + u 

+ V3 + 4k|2km(3 + k)(k(2m- 9) - 9) + 9(1 + re) (3 + 4k) (3 + ku 2 ) In 7~ j^j • (TV.26) 

The boundary action is given by 

/ 2 \ 3 ttTV 3 /• / f JJ \ 

5 boun dary = - *r-± Hm / (it dx J - a' x d x + ■ ■ ■) . (IV.27) 



9 «->o J \ u 

Substituting Eq. PV.26P into the boundary action, we get the correlator, the conductivity and the 
diffusion constant: 

6^ = 0^*1^- ^+^ , (I V.2 8 ) 
3 6 VT+^ 



= ,^: n - / „ 2 ; , . (iv.so) 



3 (l + re/3) 3 (l-re/3) 
8~T (1 + k)(1 + 4k-k 2 ) 

Thus, A is finite and Dr — on the critical line re = 3. 



V. DISCUSSION 



A. A remark on model H 



We have shown that single R-charge black holes belong to model B. Model B corresponds to a 
system with a conserved charge. But if there is another conserved quantity in the system, and if they 
are coupled to each other, the dynamic universality class becomes model H. For our systems, there is 
an obvious candidate for such a conserved quantity: T M „. In this case, A [defined in Eq. piI15[) ] and 
fj become singular as well, and the model is characterized by two dynamic exponents x\ and x^. 

D R = Xx^-^Xt 1 > (v.i) 
? 7 = €*> . (V.2) 

These dynamic critical exponents satisfy the following constraint: 7 

x x + x v = 4 - d - Tj . (V.3) 

However, we found that A remains finite for single R-charge black holeSjWhich suggests model B. 
Also, it has been shown that fj remains finite for the D = 5 black hole [ToL fill] . 

This is because we focus on linear perturbations. The difference between model B and H comes 
from a "mode-mode coupling," which arises from nonlinear terms in hydrodynamic equations. As far 
as linear perturbations are concerned, there is no distinction between model B and H. 



7 Equation IIV.3I I corrects the sign in front of r\ in Ref. [||. 
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B. The value of the dynamic exponent z 

Model B satisfies the dynamic scaling relation z = 4 — 77 from the theory of dynamic critical 
phenomena. What is the value of the exponent for single R-charge black holes? Let us focus on the 
D = 5 case. The value of z requires the knowledge of r\. As discussed in App. [A] r\ = 1 from static 
scaling laws, which implies z — 3. 

There is a caveat though. We extracted the exponents (i>, n) from a static hyperscaling relation, but 
the hyperscaling relation is often trustable below the upper critical dimension. Substituting the values 
of static exponents Eq. (|III.14|) into the Ginzburg criterion (|II.13|) . one gets d c = 3. So, our system is 
marginal, and the hyperscaling relation may not be trustable. One should not take the value z = 3 
too literally. 

In order to avoid this problem, it is desirable to obtain the relaxation time in the critical regime, 
which is expected to behave as 1 /r q cx q z . In the gravity side, this amounts to solve the scalar mode 
of the bulk U(l) field, and the relaxation time is given by the lowest quasinormal frequency wqnm(<7) 
of the scalar mode: 

1/t- 9 ~-3[^ QN m(<z) ] ■ (V.4) 

Thus, the dispersion relation of the quasinormal frequency should behave as 

3 [ wqnm(?) ] oc q z . (V.5) 

It would be interesting to check if the scalar mode computation really gives z — 3 implied from the 
static scaling relations [HI]. Note that only even powers of q appear in the equations of motion, so 
one expects that only even powers of q also appear in the dispersion relation as well unless analyticity 
is broken. 



C. Relation to QCD 

Our systems are single R-charge black holes, which are not realistic from the point of view of QCD. 
Moreover, according to lattice results, the finite temperature transition from hadronic phase to QGP 
(at zero baryon chemical potential) is not a phase transition but rather a smooth crossover. Thus, one 
may wonder if the type of phenomena considered in this paper has any relevance to the real QCD. 

However, the finite density transition at T = is believed to be a first-order phase transition, so 
there must be an end point of the first-order phase transition line somewhere in the phase diagram, 
which is a critical point. The precise location of the critical point is still unknown, and one of major 
goals of future RHIC experiments is this critical end point search (See, e.g., Ref. [lal). Physics 
near the critical point has been studied in Refs. [13, [3, [H, H3] based on the singular behavior of 
thermodynamic functions. The dynamic critical phenomena have been studied as well [U, [2^, HH ; it 
was argued that QCD belongs to model H with the dynamic critical exponent for the baryon diffusion 
constant z ~ 3. 



D. Future directions 

We have analyzed the dynamic critical phenomena for single R-charge black holes in the linear 
regime. One interesting direction is to analyze the system in the nonlinear regime. In this case, the 
theory of dynamic critical phenomena predicts (genuine) model H. It would be interesting to see if the 
black hole computation is consistent with model H predictions. A particularly interesting quantity 
is the exponent x n in Eq. (|V.2[) ; if x n > 0, the shear viscosity diverges as well. On the other hand, 
according to the AdS/CFT duality, fj/s = l/(47r) in the strong coupling limit for all known examples. 
Thus, x v > means the violation of the universality for fj/s. (However, it still satisfies the conjectured 
universal bound for fj/s [24| . and it is known that fj has only a weak critical singularity.) 
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In an unstable system, the ordering process occurs and the system tends to change to an ordered 
phase (spinodal decomposition). Because of the second-order phase transition, a single R-charge 
black hole solution should be replaced by a new bulk solution which corresponds to a new phase. It is 
important to find such a bulk solution and to analyze the dynamics after the transition. A particularly 
important question is the dynamics of the interface between two domains. When a system with second- 
order phase transition is quenched, phase separation occurs and domains are formed with the interface. 
The interface is described by a kink solution which connects two phases. It would be interesting to 
find and to analyze the bulk solution which corresponds to the kink. 

In this paper, we have considered only single R-charge black holes, but it would be interesting to 
extend the analysis to the other black holes. Because there is ample evidence that fj/s = l/(4ir) in 
the linear regime, it would be better to consider the other transport coefficients such as the charge 
diffusion constant, heat conductivity, and the bulk viscosity. Also, we consider only model B and H, 
but it would be interesting to identify black holes with the other dynamic universality classes. 

The analysis such as the one given in this paper may be useful to condensed-matter applications of 
the AdS/CFT duality. Namely, given a condensed-matter system, the static and dynamic universality 
classes may be useful to identify the dual black holes. For example, it would be interesting to find a 
bulk black hole solution in the same universality class as the Ginzburg-Landau theory [i.e., the same 
set of critical exponents as Eq. (|II.12p ] . 



We would like to thank Hirotsugu Fujii, Kenji Fukushima, Teiji Kunihiro, Chiho Nonaka, and Hi- 
rofumi Wada for useful discussions. This research was supported in part by the Grant-in-Aid for 
Scientific Research (20540285) from the Ministry of Education, Culture, Sports, Science and Technol- 
ogy, Japan. 

APPENDIX A: STATIC CRITICAL EXPONENTS FOR THE D = 5 SINGLE R-CHARGE 

BLACK HOLE 

The static critical exponents for the single R-charge black holes were obtained in Refs. [8, 9], but 
we include the D = 5 computation here for completeness. 

First, let us approach to the critical line (T c , jjL c ) with /i = fi c . Denote the deviation from the critical 
line by tr- Since 7rT = /i c \/(l + k)/2k, 
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-2 cx (2 - k) 2 . 



(A.l) 



Then, thermodynamic quantities behave as 




N 2 cl xl (1 + s)(2 + s)(3-k) 1 



-1/2 



(A.2) 



P - Pc 



1/2 



(A.3) 



XT 



16 7T 2 k(2-k) 2-k t 



-1/2 



(A.4) 



which determine 3 static critical exponents 



a = (3 



7 = 1/2. 



(A.5) 



1G 



Similarly, let us approach the critical line with T = T c , and denote the deviation from the critical 
line by Since T = 2T C \/T+k;/(2 + k), 

e„:=l---.= .oc(2- K ) 2 . (A.6) 

Mc ( 2 + K )(V^+V2) 



The R-charge density then behaves as 

p _ pc _ ^T>Ml + ^-H^)> , .p.,, , _ e y, , (A . ?) 

which determines 5 as 

(5 = 2. (A.8) 

The static critical exponents (|A.5p and (|A.8|) satisfy the static scaling relations 

a + 2/3 + 7 = 2, 7 = 0(5-1). (A.9) 

The static exponents v and 77 require the R-charge correlator, but it has not been evaluated. Instead, 
if we use the scaling relations 

2-a = vd, 7 = 1/(2-77), (A.10) 

we get 

M=(l, 6 -^), (A.11) 



or for d = 3, 



2d ' 3 



(/'.//) = (i,l) ■ (A. 12) 



Note that the value of 77 differs from the usual mean-field theory value 77 = 0. 
Consequently, we get 



in. /., (1,1,1,2,1,1) . (A. US) 
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